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We analyze the laser cooling of polarizable particles by continuous dispersive position detection 
and active feedback. The magnitude of the dissipative force is proportional to the particles' photon 
scattering rate into the detector, while its velocity dependence is determined by the programmable 
frequency dependence of the loop gain. The method combines final temperatures near the recoil 
limit with large velocity capture range, and is applicable to multilevel atoms or molecules. 



PACS numbers: 32.80.Pj 

Low-temperature, high-brightness sources of particles 
are of interest in many areas of physics. In particle ac- 
celerators, stochastic cooling has been a key technology 
to increase collision rates lj. In a single-particle pic- 
ture, an electrical signal that depends on a particle's 
motion generated in one part of the accelerator ring is 
used to apply a correction force further downstream. In 
view of Liouville's theorem, which forbids compression 
of phase space volume by conservative forces, the cool- 
ing of a finite-number sample has been explained as an 
expulsion of the empty phase space between occupied 
points pj. Feedback cooling has been applied to other 
systems including a resonant gravity gradiometer ]2j , a 
single mirror vibration mode via radiation pressure lij-Llf ■ 
and a single axially oscillating electron in a Penning trap 
|{| Feedback concepts have been further extended to 
the quantum domain 0, IS IS 13 • I n or der to increase 
the cooling bandwidth, an optical version of stochastic 
cooling has been proposed for accelerators [ljj ■ 

For atomic samples, Raizen et al. have proposed 
to optically measure fluctuations of the center-of-mass 
momentum, and to apply momentum kicks with laser- 
induced dipole forces [lj. Balykin and Letokhov have 
suggested a cooling scheme based on velocity measure- 
ment and separation for individual atoms |T^ . Recently, 
Steck et al. have analyzed quantum feedback cooling of 
a single atom strongly coupled to a resonator [l3|. In 
cavity experiments , si ngle atoms have been trapped by 
applying feedback [lj, [l5|] , while in an optical lattice, a 
sample's center-of-mass oscillation has been damped [l6j| . 

In this Letter, we analyze the laser cooling of a gas 
of polarizable particles by continuous dispersive position 
measurement and optical feedback. In particular, the 
sign and velocity dependence of the dissipative force are 
determined exclusively by the sign and frequency depen- 
dence of the loop gain. The force is proportional to the 
radiation pressure associated with photon scattering into 
the detector, but is otherwise independent of the target 
level structure. We derive closed-loop heating rates due 
to light and photodetector quantum noise, and due to the 
thermal motion of the sample. Using an optical resonator 
for signal enhancement in combination with nearly shot- 
noise-limited photodetection, a large variety of atomic or 
molecular samples can be cooled to temperatures near 



the photon recoil energy. 

In the following we adopt the single-particle viewpoint 
of stochastic cooling lj, where the signal generated by 
a single atom is propagated through the linear feedback 
loop, and the thermal motion of the other atoms in the 
sample constitutes a source of heating. Consider an 
atom of mass m moving as x a =vt at constant velocity v 
in a weak periodic potential V(x,t)=U(t) cos2fca; whose 
depth U(t) <C mv 2 /2 is influenced by the atom's mo- 
tion. Since the time variation of the force on the atom 
f(t)=2kU(t) sm2kvt arises both from the spatial varia- 
tion cos 2kx of the potential, and from the dependence of 
the depth U(t) on the atom's trajectory, a frequency com- 
ponent of U(t) in phase with the atomic-motion-induced 
force variation sin 2kvt produces a non-zero average force 
that can heat or cool the atom. To establish this type 
of dissipative system, we propose to use a standing light 
wave to both monitor the position x a (t) of the polariz- 
able atom by means of its index of refraction [lj, [lTJ > 
and to generate a periodic potential of adjustable depth 
U(t) via the light shift. 

Viewed in the time domain, the atom is cooled because 
it finds itself moving uphill on a deeper potential than ex- 
perienced when moving downhill (see also Refs. (l7tll8| L 
Of course, this is only true for an appropriately chosen 
transfer function mapping x a (t) onto U(t). In the fre- 
quency domain picture, both the moving atom and the 
feedback actuator modulate the light wave at the Dopplcr 
frequency 2kv. Under appropriate conditions, amplitude 
and phase modulation from both processes interfere to 
produce asymmetric sidebands in the light exiting the 
system. If the blue Doppler sideband is stronger, the 
atom will be cooled at the rate at which energy is car- 
ried away by the frequency-shifted light. The atom- light 
interaction, and consequently the cooling force, can be 
increased by means of an optical resonator (Fig. 

The resonator of finesse F and field decay rate constant 
7 C supports a standing-wave Gaussian TEMqq mode 
of wave number k with waist w. Incident light pro- 
duces an electric field of amplitude 2E C at an antin- 
ode (x=0) inside the resonator. An atom with com- 
plex polarizability a moving on the resonator axis co- 
herently scatters photons into free space at an average 
rate r sc =fc 3 |i?e(a:)-E c | 2 /(67reo?i), and experiences an op- 
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FIG. 1: Schematic of continuous-feedback cooling. The mov- 
ing atom's effective index of refraction periodically detunes 
the standing- wave resonator by 5 at (t) - The fractional change 
in cavity transmission e(t) acts as an error signal that is ap- 
plied to change the input power by a fraction £ib(t). The 
feedback open-loop gain H(s) is set by an external electronic 
circuit. The input light is initially detuned by Si from cavity 
resonance. 



tical potential E/oCOs2fca; of depth Uq—— \E c \ 2 Re(a)/2. 
Conversely, the forward-scattered radiation by the spa- 
tially varying induced atomic dipole, or equivalently, 
the atom's index of refraction, results in an atomic- 
position-dependent detuning <5 a t=C7c cos2fcx a of the res- 
onator from its average resonance frequency u> c |l7l Il8| . 
Here (^=hr]T sc /Uo is a dimensionless parameter char- 
acterizing the atom-cavity coupling, and the quantity 
rj=6F/ (irk 2 w 2 ) can be interpreted as the fraction of pho- 
tons scattered into one direction of the resonant cavity. 

For incident light detuned by an amount Si=ck—u> c rel- 
ative to the resonator, the time- varying detuning <5 at (i) of 
the resonator by the moving atom changes the resonator 
transmission by a fraction e(t). The transmitted power is 
measured and used as an error signal in a feedback loop 
to adjust the incident power by a fraction — £fb(t} (Fig. 

. In order to ignore cavity- induced forces 0, , we 
take the resonator linewidth 2j c to be much larger than 
both the feedback bandwidth and the Doppler frequency 
2kv, such that the intracavity power P c (t) adjusts in- 
stantaneously to a value determined by the total light- 
resonator detuning St(t)=Si — S a t(t). Then the fractional 
deviation e(t)=P c (t) /P — 1 of the intracavity power from 
its unperturbed value Pq (for <5 a t=£fb=0) is given by 



s(t) 



7 C 2 + ^ 



(l-£fb(0) - 1- 



(1) 



Assuming |<$ a t|/7c £fb <C 1, we can linearize the fractional 
change e{t) in resonator transmission, 



e{t)=r 1 - 1 5 &t (t)~e ih {t). 



(2) 



In this approximation, the moving atom modulates the 
intracavity power by an amount rS^t/^/c proportional to 
the normalized resonator slope r=25[-f c / (j 2 + 8 2 ), while 
the feedback loop adjusts the incident power by a frac- 
tion — £fb(i). If the atom's kinetic energy far exceeds 
the light shift Uq, then to lowest order 5 at (t) is deter- 
mined by the atom's unperturbed motion x a —vt. We in- 
troduce the open-loop feedback gain H(s) for the Laplace 
transformed quantities ifb(s),i(s) via S{b=He. Then 
H(iu))=H\(uj) + iHiiuS), with real and imaginary parts 



H\ and H 2 , respectively, is the complex gain in the fre- 
quency domain. If the loop is stable, the steady-state 
solution in the time domain is given by 



e(t) - r(- 



'l + Hi{2kv)) cos2kvt + H 2 {2kv) sin2kvt 
\1 + H(2ikv)\ 2 



(3) 

Hi(2kv) and H 2 (2kv) are the open- loop gain in phase 
and in quadrature with the atomic-motion-induced inten- 
sity modulation <5 a t(i)=C7c cos 2kvt, respectively. They 
determine the closed-loop signal e(t), and thereby the 
time variation of the optical-potential depth U(t)=Uo(l+ 
s(t)). In the limit Uq <C mv 2 /2, work is done on the 
atom at a rate W=e(t)f u (t)v to lowest order, where 
/ u = 2kUo sin 2kvt is the unperturbed force. The com- 
ponent sin2/ci>£ of e(t) in phase with / u produces a dis- 
sipative force, whose spatial average f=(W)/v can be 
written as 



f(v) = hkrjTs 



rH 2 {2kv) 
\1 + H{2ikv)\ 2 



(4) 



This expression, derived without invoking the two-level 
or rotating-wave approximations, is valid for arbitrary 
laser detuning from atomic resonances below saturation. 
The friction force / depends on target parameters exclu- 
sively through the scattering rate T sc and is controlled 
by the average intracavity intensity. / is proportional 
to the rate of photon momentum transfer 2hkr]T sc due 
to backward scattering into the resonator at rate 77F SC , 
multiplied by a dimensionless function of the atomic ve- 
locity. In particular, the sign and velocity dependence 
of / are determined by the frequency-dependent loop 
gain H{iu>). They are independent of the atom's level 
structure, such that atoms in different internal states, or 
different species, can be cooled simultaneously. Cooling 
occurs if the product of the resonator line slope r and the 
quadrature loop gain H 2 {2kv) at the Doppler frequency 
2kv is negative. The force is maximized at a resonator- 
light detuning <5i=±7 c that gives the largest slope r=±l. 
In the following we assume <5j=— "f c (r=— 1), such that 
H(s) > corresponds to negative feedback. 

The cooling force / is proportional to the quadrature 
component sin 2kvt of the intracavity light modulation in 
closed loop, given by H 2 /\l + H\ 2 . In particular, for very 
small or very large open-loop gain (\H\ 2 =H 2 +H 2 -C 1 or 
I if 1 2 3> 1), this relevant quadrature of the intensity vari- 
ation, and hence /, will be small. The velocity-dependent 
term takes on its maximum value 1/(2 + 2Hi) when the 
quadrature gain H 2 and the in-phase gain Hi are related 
by |if 2 |=|l + -ffi|- When < 1, the feedback loop 

regeneratively amplifies the intensity variation caused by 
the moving atom, so that the force is increased by a factor 
|1 + -Hil -1 S> 1 compared to the case of small in-phase 
gain I Hi I -C 1. However, the heating of the atom due 
to noise amplification is then also increased, as we show 
below. 

The simplest stable cooling loop is a differentiator , 
Hd(iw)=iuj / (2ku) , with unity gain frequency 2ku. (The 
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FIG. 2: Feedback cooling force versus atomic velocity for 
various loop gains. The feedback loop gain can be chosen 
(as in a, b, or c) to produce a much larger velocity capture 
range than conventional Doppler cooling (d) , while still yield- 
ing the same friction coefficient near zero velocity (see inset). 
The open- loop gains plotted are (a) H(s)=s (differentiator), 
(b) ff(s)=s(l + s/10)/(l + 8s/10), ( c ) H(s)=s(l + s/W)(l + 
s/100)/((l + 8s/10)(l + s/20)), and (d) H(s) = s(l + s/2) 
reproducing the Doppler force. Here s=iv/u, where u sets 
the unity gain frequency for the differentiator loop. 



gain can be rolled off at high frequency outside the veloc- 
ity range of interest to the cooling.) We further assume 
for simplicity that T ac is made independent of atomic ve- 
locity by detuning the laser frequency from atomic res- 
onances by an amount much greater than the Doppler 
shift 2kv. The differentiator-induced cooling force /d is 
then given by 



fd(v) = -hkr]T s 



(5) 



The maximum force hkrjT sc , attained for the unity- 
gain velocity v=u, is the same as in Doppler cooling 
at the photon scattering rate ?7r sc , but the differentia- 
tor force, falling off as 1/v rather than 1/w 3 , provides 
a substantially larger velocity capture range than the 
Doppler force. Fig. [2] also shows that higher-order loops 
can substantially further increase the velocity capture 
range while maintaining the low-velocity friction coeffi- 
cient ao=df /dv that determines the final temperature. 
«o can be increased in real time simply by increasing the 
loop gain H, i.e. by reducing u. We also note that a loop 
H(s)=(s/T')(l+(s/(2T')) reproduces the velocity depen- 
dence of the Doppler force |l9j for a (fictional) transition 
with linewidth 2T' < 2j c . 

Technical or quantum noise in the light-induced po- 
tential U(t) will heat an atom moving at velocity v in 
proportion to the spectral noise density at the mod- 
ulation frequency 2kv of the unperturbed force. We 
define the single-sided fractional spectral noise den- 
sity S c (ui)—(2/tt) Jp 00 dt(e(0)e(t)} cosujt, normalized such 
that J °° duiS c (ui)=(e 2 (0)} is the mean-square fractional 
intensity noise, and calculate the heating power due to 



momentum diffusion [Tsj . 



Wfluct 



nk 2 



U 2 S c {2kv). 



(G) 



This general formula describes the heating of an atom 
moving at velocity v due to classical intracavity intensity 
noise. For quantum noise, Eq. needs to be multiplied 
by a factor of two to take into account fluctuations of the 
induced atomic dipole |l9j |. Assuming negligible techni- 
cal noise, the shot noise of the light incident on the cavity 
with power P; corresponds to a fractional power spectral 
density given by Sf sn =hck/(Pin). In closed loop, the in- 
tracavity noise density is S , c psn (w)=5 , i psn /|l + H(iuj)\ 2 for 
u> <C 7c. In addition, shot noise in the detector pho- 
tocurrent at limited quantum efficiency q < 1 arising 
from the undetected photons will produce a false error 
signal, causing the feedback loop to generate (classical) 
intracavity intensity noise. This photoelectron detection 
noise corresponds to an uncorrelated fractional noise den- 
sity of magnitude S' osn =S' i psn (<7~ 1 — 1) in open loop, and 
S^=S esn \H\ 2 /\l + H\ 2 in closed loop. The total shot 
noise from these two sources S^ n =S^ sn + S° sn , expressed 
in terms of the intracavity power P c =PiF j '(2ir) for the 
detuned cavity with |J;|=7 C , is given by 



Fhck i + lif^por 1 - 1) 



2^ 2 P, 



(7) 



Eq. JfJJl then yields a closed-loop heating rate W sn due 
to shot noise 



W S1 



2+\H(2ikv)\ 2 (q- 1 - 1) 
: \1 + H(2ikv)\ 2 : 



(8) 



where E r —h 2 k 2 /(2m) is the recoil energy. We see that 
Wsn can be expressed as the recoil heating rate 2i? r ?7r sc , 
multiplied by a dimensionless function that depends on 
the loop gain H and the detector quantum efficiency q. 
In the absence of feedback (H=0), Eq. (JHJ) reproduces 
the well-known recoil heating due to scattering into 
the resonator, where each backscattering event, occur- 
ing at a rate (?7/2)r sc for the detuned cavity, heats the 
atom by an amount AE r . For unity quantum efficiency 
(g=l) the feedback loop reduces (increases) this heating 
by suppressing (enhancing) the intracavity fluctuations 
for |1 + H\ > 1 (|1 +H\< 1). For q < 1, the feedback 
loop also generates intensity noise originating from ran- 
dom photoelectron detection. In addition to the cavity 
heating W sn , scatte ring into free-space heats the atom at 
a rate W is =2E r T sc 

For the differentiator loop H^iito) — iuj/(2ku), a unity- 
gain velocity u=(q~ 1 — l + 2r/~ 1 )hk/m minimizes the tem- 
perature Td, yielding 



k B T d = AE r 



(l+^jq- 1 -1 + 2TT 1 ) 

r) 



(9) 



where fcp is Boltzmann's constant. Thus for good atom- 
cavity coupling, r) = 1, the final temperature is limited 
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by the finite detection efficiency to 8(l+q~ 1 )E r /kB- This 
implies that except for the lightest atoms and molecules, 
microkelvin temperatures can be reached even for limited 
detector quantum efficiency q. 

For a thermal sample consisting of N+l atoms at tem- 
perature T, any chosen probe atom will be heated by 
the intracavity intensity noise induced by the other N 
randomly moving atoms. The corresponding closed-loop 
spectral noise density S^(u>), expressed in terms of the 
sample's thermal velocity Vt^^sT/m) 1 ' 2 , is 



N( 2 exp(-uj 2 /(8k 2 v? h )) 
8^kv th \l + H(ioj)\ 2 



(10) 



resulting in a collective heating rate Wn for an atom with 
velocity v given by 



Wk(v) 



E rV T s 



2-KNrir s 



\1 + H{2ikv)\ 2 2kv t h 



• exp 



As is typical of stochastic cooling, the heating over- 
whelms the cooling at too large atom number N or too 
large cooling speed, i.e. photon scattering rate T sc . For 
a differentiator loop, we find that the net cooling power 
averaged over the thermal sample, (W e j)=(/d' y + Wn) 
is maximized for u w t>th- The cooling rate constant, 
7d=— 2(Wd) / {k-eT) can then be written in the form 



7d 



fc^th 
67V 



(2r cav r cav ) 



(12) 



where f ca:v =2Nr)T sc h/ (k^T) is the total scattering rate 
into the cavity normalized to the sample temperature. 
The cooling rate is maximized for r cav =l, yielding a 
rate constant 7d = kv^/i&N). We see that the thermal 
Doppler broadening kvth takes the role of the stochastic- 
cooling bandwidth 1J, and that the cooling of smaller 
subsamples can proceed faster 0, |n| . For a sample of 
iV=10 8 magnetically trapped CaH molecules at 0.4 K 
pof cooled with light of wavelength 760 ran, the opti- 
mum cooling rate is 7d=0.1 s _1 , attained at a photon 
scattering rate ?7r sc =3 x 10 2 s _1 . (Here we are assuming 
that the different degrees of freedom are mixed 0, e.g., 
by the trapping potential.) An ensemble of 10 6 molecules 



at room temperature can be cooled at a rate 7d=340 s _1 
for a photon scattering rate of r;r sc =3 x 10 7 s _1 . Eq. IT21 
implies that the differentiator-induced cooling of a ther- 
mal beam containing N atoms with light of wavelength 
A requires a characteristic length L=6NX/tt, indepen- 
dent of the initial velocity of the beam, since the cooling 
bandwidth kv t h is larger for a faster beam. 

It may seem surprising that an error signal propor- 
tional to the atom's index of refraction generates a cool- 
ing force that is proportional to the photon scattering 
rate into the detector. The reason is that the index of 
refraction arises from the superposition of the forward- 
scattered field with the much larger incident field. The 
dispersive signal can thus be viewed as a homodyne 
method to measure the small forward-scattered field. Ho- 
modyne detection, however, cannot be used to improve 
over the shot-noise limit resulting from direct detection 
of the scattered photons. The role of the resonator is to 
enhance the detection signal-to- noise ratio [12| , or equiv- 
alently, the photon scattering rate ?7r sc into the detector. 
In the absence of the resonator, rj is replaced by the detec- 
tion solid angle Af2=3/ (k 2 w 2 ), resulting according to Eq. 
@ in much higher temperatures. This may explain why 
stochastic cooling has not been observed in free-space ex- 
periments using laser beams with large waist w 3> fc _1 
for detection [llL[T^|. 

In conclusion, we have derived simple analytical ex- 
pressions for the cooling and heating of a gas interact- 
ing with a laser beam whose intensity is adjusted in re- 
sponse to the atoms' motion. If a resonator is used for 
signal enhancement in combination with photodetection 
near the quantum limit, microkelvin temperatures can 
be reached. By appropriately designing the frequency- 
dependent loop gain, large friction coefficients can be 
combined with a very large velocity capture range. The 
friction force scales with the photon scattering rate by the 
target atoms, but is otherwise independent of the target 
level structure. Therefore the proposed method repre- 
sents a promising route for laser-cooling new atomic or 
molecular species. 
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